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ABSTRACT 

Aims. We use advanced 3D NLTE radiative magnetohydrodynamic simulations of the solar atmosphere to carry out detailed tests of 
chromospheric diagnostics at millimeter and submillimeter wavelengths. 

Methods. We focused on the diagnostics of the thermal structure of the chromosphere in the wavelength bands from 0.4 mm up to 
9.6 mm that can be accessed with the Atacama Large Millimeter/Submillimeter Array (ALMA) and investigated how these diagnostics 
are affected by the instrumental resolution. 

Results. We find that the formation height range of the millimeter radiation depends on the location in the simulation domain and is 
related to the underlying magnetic structure. Nonetheless, the brightness temperature is a reasonable measure of the gas temperature 
at the effective formation height at a given location on the solar surface. There is considerable scatter in this relationship, but this is 
significantly reduced when very weak magnetic fields are avoided. Our results indicate that although instrumental smearing reduces 
the correlation between brightness and temperature, millimeter brightness can still be used to reliably diagnose electron temperature 
up to a resolution of 1". If the resolution is more degraded, then the value of the diagnostic diminishes rapidly. 

Conclusions. We conclude that millimeter brightness can image the chromospheric thermal structure at the height at which the 
radiation is formed. Thus multiwavelength observations with ALMA with a narrow step in wavelength should provide sufficient 
information for a tomographic imaging of the chromosphere. 
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1. Introduction 

As a result of its particularly complicated physics, the chro¬ 
mosphere remains the least understood layer of the solar atmo¬ 
sphere. The past decade of modeling efforts has enormously ad¬ 
vanced our understanding of this enigmatic layer, however. In 
recent years, 3D modeling of the chromosphere, including much 
of the necessary physics, has finally started to become feasible. 
Advanced 3D magneto-hydrodynamic models run from the up¬ 
per convection zone through the chromosphere into the lower 
corona (Gudiksen et al. 2011) and take into account physical ef¬ 
fects such as magnetic field dynamics, thermal conduction, am- 
bipolar diffusion, and the Hall effect, non-equilibrium ionization, 
and NLTE effects (see, e.g., Carlsson 2009, for a review). These 
simulations provide the opportunity to model the formation of 
various chromospheric lines and continua and to study the diag¬ 
nostic potential of chromospheric data (Leenaarts et al. 2013a,b; 
Pereira et al. 2013; Stepan et al. 2012; Leenaarts et al. 2012). 


ter and millimeter wavelengths (Loukitcheva et al. 2004). The 
intensities of submm/mm continua, whose source function can 
be treated in LTE (Rutten 2003), depend linearly on temperature 
and therefore may be able to provide a sensitive test of numeri¬ 
cal models. The additional advantage of submm/mm continua is 
that their intensities can be easily synthesized from the models. 

Loukitcheva et al. (2004) compared a large collection of sub¬ 
millimeter and millimeter observations with synthetic brightness 
temperatures calculated in classical standard models by Fontenla 
et al. (1993) and the ID dynamic model of Carlsson & Stein 
(1995) and demonstrated that both the classical and the dynamic 
model provide a reasonable fit to observed temporally and spa¬ 
tially averaged millimeter data. The analysis of the ID dynamic 
simulations of Carlsson and Stein further revealed that radio 
emission at millimeter wavelengths is extremely sensitive to dy¬ 
namic processes in the chromosphere if these are spatially and 
temporally resolved (Loukitcheva et al. 2004, 2006). 

Wedemeyer-Bohm et al. (2007) were the first to study 
submm/mm brightness synthesized from 3D radiation (magneto- 
hydrodynamic simulations, using those of Wedemeyer et al. 
(2004) and Leenaarts & Wedemeyer-Bohm (2006). These 3D 
(M)HD models were constructed under the assumptions of LTE 
for the radiative transfer and the equation of state to make 
the computations feasible. In addition, Wedemeyer-Bohm et al. 


Chromospheric diagnostics, such as Ca n H&K, H a, Ca n 
854.2 nm, and Mg n h&k lines, suffer from the fact that they 
are formed out of local thermodynamic equilibrium (LTE) and 
thus decouple from the local conditions. In addition, their for¬ 
mation is complex. An alternative chromospheric diagnostic is 
provided by observations of the radio continuum at submillime- 
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Fig. 1 . Horizontal cuts through a 3D MHD simulation snapshot. Plotted are the electron temperature (top panel), electron number density (middle 
panel), and magnetic field strength (bottom) at geometrical heights of 650, 900, 1500, 1700, and 2000 km. To enhance the temperature contrast, 
the displayed temperature range is set to that of the synthetic brightness maps shown in Fig. 2. The range of temperatures is given in the lower 
right corner of the frames. The field size is 24 Mm x 24 Mm. 


(2007) studied submm/mm brightness from representative model 
snapshots with non-equilibrium electron densities and a weak 
magnetic field. The authors reported mm brightness maps with 
filamentary brightenings, resulting from shock-induced thermal 
structure, and fainter regions in between. The average brightness 
temperature and relative contrast were found to increase with 
wavelength. 

This paper is the first in a series in which we use advanced 3D 
radiative magnetohydrodynamic simulations of the solar atmo¬ 
sphere performed with the Bifrost code (Gudiksen et al. 2011) 
to carry out detailed tests of chromospheric diagnostics at mil¬ 
limeter and submillimeter wavelengths. In this paper we focus 
on the diagnostics of the chromospheric thermal structure in the 
wavelength bands accessible to the advanced radio interferom¬ 
eter Atacama Large Millimeter/Submillimeter Array (ALMA). 
Located on the Chilean Chajnantor plateau above 5000 m alti¬ 
tude, ALMA consists of a giant array of 12 m antennas (the 12 
m array) with baselines of up to 16 km, which aims at achiev¬ 
ing high spatial resolution, and an additional compact array of 
7 m and 12 m antennas that greatly enhance ALMA’s ability to 
image extended targets. At the time of writing, ALMA observes 
at wavelengths in the range 3.6 mm to 0.4 mm (84 to 720 GHz) 
in dual polarization mode with prospects of going to longer (up 
to 10 mm) and shorter (down to 0.3 mm) wavelengths in the 
future. We simulate the mm/submm brightness at ALMA wave¬ 
lengths from the 3D MHD simulations and compare the results 
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with the actual temperature distributions to assess the accuracy 
of the mm-wavelength diagnostics. 

The structure of the paper is as follows: In Sects. 2 and 3 we 
describe the 3D model atmosphere based on the Bifrost code and 
the submm/mm radiative transfer computations. The analysis of 
the synthetic brightness is presented and discussed in Sect. 4. In 
this section we also investigate the effect of spatial smearing of 
the model brightness in the light of future interferometric obser¬ 
vations with ALMA. We summarize our findings and conclude 
in Sect. 5. 


2. 3D Model atmosphere 

We studied the formation of millimeter and submillimeter con- 
tinua in a snapshot of a 3D radiation-MHD simulation performed 
with the Bifrost code (Gudiksen et al. 2011). We used snapshot 
385 of the simulation “en024048_hion” made available through 
the Interface Region Imaging Spectrograph (IRIS, De Pontieu 
et al. 2014) project at the Hinode Science Data Centre Europe 
(http://sdc.uio.no). The same snapshot was used by Leenaarts et 
al. (2012) to investigate the H a line formation and by de la Cruz 
Rodriguez et al. (2013) to study Ca n 8542 spectra. It was also 
used by Leenaarts et al. (2013a) and other papers in that se¬ 
ries to investigate the formation of the Mg ii h and k lines in 
the solar atmosphere. Bifrost solves the equations of resistive 
MHD on a staggered Cartesian grid. The simulation includes 
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optically thick radiative transfer in the photosphere and low 
chromosphere, parameterized radiative losses in the upper chro¬ 
mosphere, transition region, and corona (Carlsson & Leenaarts 
2012), thermal conduction along magnetic field lines (Gudiksen 
et al. 2011), and an equation of state that includes the effects of 
non-equilibrium ionization of hydrogen (Leenaarts et al. 2007). 

The simulation covers a physical extent of 24 x 24 x 16.8 
Mm, with a grid of 504 x 504 x 496 cells, extending from 
2.4 Mm below the average height of T 500 = 1, which corre¬ 
sponds to Z = 0, to 14.4 Mm above, so that it covers the up¬ 
per convection zone, photosphere, chromosphere, and the lower 
corona. The horizontal axes have an equidistant grid spacing of 
48 km (0.06"), the vertical grid spacing is nonuniform, with a 
spacing of 19 km between Z = -1 Mm and Z = 5 Mm. The 
spacing increases toward the bottom and top of the computa¬ 
tional domain to a maximum of 98 km. The simulation contains 
a magnetic field with an average unsigned strength of 50 G in the 
photosphere, concentrated in the photosphere in two clusters of 
opposite polarity lying 8 Mm apart, representing two patches of 
quiet-Sun network. We refer to the snapshot from this simulation 
we analyzed simply as the 3D MHD snapshot. For more details 
of this snapshot we refer to Carlsson et al. (2013, 2015). 

In Fig. 1 we show examples of horizontal slices through the 
3D MHD snapshot of electron temperature (upper panel), elec¬ 
tron number density (middle panel), and magnetic field strength 
(bottom) taken at the geometrical heights Z = 650, 900, 1500, 
1700, and 2000 km. These heights very roughly correspond to 
the effective formation heights of the radiation at 0.4, 1.0, 3.0, 
4.5, and 10 mm, respectively (see Sect. 4.2). Between the mag¬ 
netic patches the temperature shows filamentary structure that 
follows the general orientation of the magnetic field. A hot chro¬ 
mospheric canopy lies above and between the photospheric mag¬ 
netic field concentrations, and coronal temperatures are reached 
already at heights of around 1700 km. The electron density 
displays a similarly complex pattern, with a clear anticorrela¬ 
tion with the temperature above magnetic patches at heights 
above 1700 km. The magnetic field concentrations expand with 
height from their photospheric footpoints (with a maximum field 
strength of 2200 G) to a volume-filling field higher up with a 
maximum strength of 60 G at a geometrical height of 2000 km 
(see also Fig. 5). 

3. Synthesis of mm brightness 

We calculated the submillimeter and millimeter radiation emerg¬ 
ing from the 3D MHD snapshot at 32 wavelengths from 0.4 mm 
to 10 mm. In the calculations we assumed each vertical column 
in the 3D snapshot to be an independent ID plane-parallel atmo¬ 
sphere. The radiative transfer was computed under the assump¬ 
tion that bremsstrahlung opacity is responsible for the mm con¬ 
tinuum radiation in the quiet Sun. We distinguished two sources 
of opacity: opacity due to encounters between free electrons and 
protons, and opacity due to encounters between free electrons 
and neutral hydrogen. They are commonly referred to as H° 
opacity and H" opacity. The corresponding absorption coeffi¬ 
cients depend on the effective number of collisions undergone 
by an electron per unit time, including collisions with protons 

(Xep K and with neutral hydrogen atoms (Xcu N ‘% - )• 

Expressions for calculating the radiative transfer were taken 
from Zheleznyakov (1996) and can also be found in Loukitcheva 
et al. (2004). In the presence of a magnetic field in quasilongi¬ 
tudinal approximation, the absorption coefficients for extraordi¬ 
nary ( e ) and ordinary (o) modes are determined by the relation 


(Zlotnik 1968) 



(1) 


where x% 0 characterizes absorption (^ ep or Xqh) in the absence 
of the magnetic field, co is the observing frequency, oj b is the 
gyrofrequency, a is the angle between the magnetic field and 
the line of sight, the minus sign corresponds to the e-mode, 
and the plus sign to the o-mode. For a weak magnetic field 

(| cos a |«: 1) the quasilongitudinal approximation is valid 

for all angles between the magnetic field and the line of sight, 
except for narrow intervals of angles close to transverse propa¬ 
gation (Zheleznyakov 1996). 

We considered the radiative transfer equation in terms of 
brightness temperature T^° separately for the e- and omode, 


djT 

dr 


_ rp rpC^O 

~ 7 e _ 7 b ’ 


( 2 ) 


where T e = T e (l) is the profile of the kinetic temperature along 
the light path, r(/) = J t Xe,o(J)dl is the optical depth, / is geo¬ 
metrical distance along the light path, and^ e ,o is the correspond¬ 
ing absorption coefficient. The solution of the radiative transfer 
equation on a discrete grid can be expressed in terms of T e as 
(Hagen 1951) 


T e b ’° = Yj(1~ exp-*» A/ ') T r e exp- (3) 

r= 1 


where Xe,o stands for Xe or Xo at position r, r = 1 corresponds 
to the layer at the top of the atmosphere, Ah is the grid distance. 
The items within the sum represent the contribution of various 
layers to the emerging radiation intensity. We refer to them as the 
values of the (unnormalized) contribution function (CF) to T ^°. 
The total brightness 7b and the degree of circular polarization P 
at wavelength A are defined as 


T e + T o 
2 


Tg _ 'T'O 

p _ £b_£b_ 

Tf + T? ’ 

b b 


(4) 


Hereafter we discuss the total brightness at millimeter wave¬ 
lengths as a primary diagnostics of the chromospheric thermal 
structure, leaving the discussion of the circular polarization for 
the next paper of the series, which will be on the subject of the 
chromospheric magnetic field. 


4. Results 

4.1. Distribution of brightness temperature 

In Fig. 2 we show the results of the brightness calculations for 
0.4 mm (ALMA band 9), 1 mm (band 7), 3 mm (band 3), 4.5 mm 
(band 2) and 10 mm (band 1). The brightness temperature at 
all considered wavelengths exhibits a complex pattern of inter¬ 
mittent bright and dark regions (top row of Fig. 2) similar to 
those seen in the electron temperature and density horizontal cuts 
through the model atmosphere (see Fig. 1). In the central part of 
the field of view the most prominent features are bright elon¬ 
gated fibrils that extend outward from the network patches and 
are aligned along the chromospheric magnetic field lines. At the 
shortest wavelengths the imprint of the central loop-like mag¬ 
netically formed structures is weak, but it becomes more pro¬ 
nounced toward longer wavelengths. 
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Fig. 2. Maps of simulated mm brightness at the resolution of the model for the wavelengths 0.4, 1.0, 3.0, 4.5, and 10 mm (labeled in the top 
row). The horizontal white line in the top row marks the position of the vertical cut at Y=0.9 Mm displayed in Fig. 5. We also plot the effective 
formation heights of the mm radiation (bottom row) and the corresponding electron temperatures, electron number densities, and magnetic field 
strength taken at the effective formation heights (three middle rows). The values of the most extreme 1% of points are clipped in the brightness 
maps. The displayed temperature range for electron temperatures is set to that of the mm brightness maps. The range of covered temperatures is 
given near the bottom of the frames. The field size is 24 Mm x 24 Mm. 


To quantify the brightness distributions at mm wavelengths, 
we list in Table 1 some statistics including the lowest and highest 
brightness temperature, average brightness temperature < 7^ >, 
RMS variation T r b ms , and the relative brightness temperature con- 

jrms 

trast r— for a number of analyzed wavelengths. The average 
brightness temperature as well as the highest brightness increase 
with wavelength as a result of the shifting contributing heights to 
higher (hotter) layers in the atmosphere (see Figs. 6, 7, 8, 9, 10, 


and 11, which are discussed in detail later). The RMS variation 
jrms a i SQ s teadily increases with A. A more complex behavior 
is displayed by the relative brightness contrast, which initially 
increases with wavelength, reaches a value of 0.22 at around 
A =3 mm, and stays almost constant at longer wavelengths. In 
Fig. 3 we plot average brightness temperatures < >, with 

T b ms as error bars, overlaid on the observed millimeter bright¬ 
ness spectra from the compilation made by Loukitcheva et al. 
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Table 1. Lowest and highest brightness temperature, average brightness 
temperature < T b >, RMS variation T r h ms , and relative brightness tern- 

jrms 

perature contrast for a number of analyzed wavelengths. 


A, mm 

jmin 

jmax 

< T b >, K 

jrms 

'j'rms 

1 b 

<T b > 

0.4 

3336 

8446 

4507 

503 

0.11 

1.0 

2706 

11094 

4786 

823 

0.17 

3.0 

2945 

13104 

6090 

1333 

0.22 

3.6 

3101 

13227 

6394 

1402 

0.22 

4.5 

3054 

13637 

6786 

1485 

0.22 

10.0 

2924 

20082 

8615 

1788 

0.21 


(2004). Gratifyingly, < Tb > calculated from the snapshot of 
the 3D MHD model (red squares in Fig. 3) are consistent within 
the error bars with the available observations of mm brightness. 
The only small discrepancy occurs at 1 mm, where, in spite of 
some overlap, the simulated 77 appears to lie somewhat below 
the observational data. Longward of 3 mm the simulated values 
are significantly higher than those reported in Wedemeyer-Bohm 
et al. (2007). 



0.1 1.0 10.0 
wavelength (mm) 


Fig. 3. Simulated millimeter brightness (red squares) overlaid on the 
observed brightness temperatures from Loukitcheva et al. (2004). Data 
obtained around solar cycle minimum, solar maximum, and at inter¬ 
mediate phases are depicted by filled circles, open circles, and stars, 
respectively. 


4.2. Effective formation heights 

The bottom row of Fig. 2 shows maps of the effective forma¬ 
tion height for five wavelengths. We define the effective forma¬ 
tion height as the height corresponding to the centroid of the 
contribution function. At 0.4 mm the radiation originates mainly 
in layers close to the traditional temperature minimum region 
and lower chromosphere with a few, very localized excursions 
to heights above 1000 km. The average formation height of this 
wavelength is around 650 km. With increasing wavelengths, we 
clearly sample higher chromospheric layers. Effective formation 
heights, averaged over all spatial locations, reach values of about 
900 km for 1 mm, 1500 km for 3 mm, 1700 km for 4.5 mm, 
and around 2000 km for 10 mm. In the expanding magnetic 
patches that overlie photospheric magnetic footpoints and show 



effective heights 



height (km) 

Fig. 4. (a) Histograms of the heights contributing more than 1% of 
the total contribution to the emerging intensity at 0.4 mm (dot-dashed), 
1 mm (dashed), 3.6 mm (solid), and 10 mm (dotted), (b) The same for 
the effective formation heights. 

enhanced temperature, we normally see down to lower chromo¬ 
spheric heights (e.g., 1200-1500 km at 3 mm), while in loop-like 
structures that connect the magnetic concentrations, we sample 
higher layers (1500-2000 km at 3 mm). Details on the forma¬ 
tion heights of individual brightness features can be found in 
Sects. 4.3 and 4.4. 

To check how representative of the true formation layers the 
effective height is, we plot in the second top row of Fig. 2 the 
electron temperature T e taken at the effective formation heights. 
For completeness, we show in the third and fourth rows of Fig. 2 
the same for the electron number density n e and the magnetic 
field strength B. Qualitatively, the properties of the atmosphere 
taken at the effective formation heights (Fig. 2) look quite similar 
to the atmospheric properties of the layers corresponding to the 
average effective heights (Fig. 1), but there is a significant differ¬ 
ence in the range of values, specifically for the magnetic field. 
From the comparison of the top row of Fig. 1, which shows T e at 
the average formation heights of mm radiation in the 3D MHD 
snapshot, with the two upper rows of Fig. 2, which depict syn¬ 
thetic Tb and T e at the effective formation heights, we conclude 
that in general, mm brightness can serve as a measure of temper¬ 
ature in the solar atmosphere, albeit not a perfect measure. We 
report the results of a correlation analysis that provides a more 
quantitative comparison in Sect. 4.5. 

Statistical distributions of the heights that contribute to the 
emerging intensity are presented in Fig. 4 in the form of his- 
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Fig. 5. XZ-profile of the 3D MHD snapshot along the cut at Y=0.9 Mm, 
marked by a horizontal line in Fig. 2. (a) Electron temperature on a log¬ 
arithmic scale clipped at 10 5 K to enhance the contrast at lower heights. 

(b) Electron number density on a logarithmic scale, (c) Magnetic field 
strength on a logarithmic scale. Overplotted dot-dashed, dashed, solid, 
and dotted lines in panels (a)-(c) depict the effective formation heights 
for the radiation at 0.4 mm, 1.0 mm, 3.6 mm, and 10 mm, respectively, 
(d) Brightness temperature at the wavelengths 0.4 mm (dot-dashed), 

I. 0 mm (dashed), 3.6 mm (solid), and 10 mm (dotted). In panels (a)- 

(c) the y-axis scale is expanded. In all panels vertical dashed lines mark 
the locations of the different structures presented in Figs. 7, 8,9,10, and 

II . 


tograms for four of the analyzed wavelengths together with the 
distributions of the effective formation heights for comparison. 
The corresponding histogram peaks in the two panels coincide 
with an accuracy of 50 km for all wavelengths except for that at 
1 mm. At 1 mm, the discrepancy between the peaks in the con¬ 
tributing and effective height histograms reaches 160 km, and 
the histograms differ in form, which suggests that statistically 
we need to be more careful at 1 mm in interpreting the effective 
formation heights. 

4.3. Analysis of a vertical slice through the snapshot 

In Fig. 5 we provide an example of an XZ-slice through the at¬ 
mosphere along the cut at Y=0.9 Mm marked in the top row of 
Fig. 2 by the horizontal white line. The panels (a-c) of Fig. 5 
display the electron temperature T e , the electron number density 
n e , and magnetic field strength B along the cut, with the effec¬ 
tive formation heights at four selected wavelengths overplotted. 
The bottom panel (d) shows the synthesized brightness along the 
cut at T=0.4 mm (dot-dashed), 1.0 mm (dashed), 3.6 mm (solid), 
and 10 mm (dotted). 

To some extent, the effective formation heights of mm radia¬ 
tion replicate the distribution of the local spatial inhomogeneities 
in electron temperature and number density along the cut. To a 
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larger extent, however, they are defined by the temperature and 
number density structure along the line of sight (in particular, the 
strong local gradients of temperature and density that influence 
opacity peaks and therefore the form of the contribution func¬ 
tion) and thus mix contributions from both high and low layers. 
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Fig. 6. Gray-scale figures of the contribution function along the cut 
at Y = 0.9 displayed in Fig. 5. Dark color corresponds to higher am¬ 
plitudes. Contribution functions in each column are normalized to the 
highest value of the column. Heights with contribution function ampli¬ 
tudes equal to 10% of the highest contribution to the emerging intensity 
are marked with blue contours. Red solid lines represent effective for¬ 
mation heights. The various panels show results for 0.4 mm (a), 1 mm 
(b), 3.6 mm (c), and 10 mm (d). In all panels the y-axis scale is ex¬ 
panded. The dashed lines are the same as in Fig. 5. 


Away from the regions with a strong magnetic field, the ef¬ 
fective formation height tends to be lower than in the vicinity 
of magnetic field concentrations (see Fig. 5c), although there are 
exceptions. Directly above the field concentrations we see a rela¬ 
tive shift of the iso-Kelvin temperature curves (and consequently 
of the effective formation heights) to higher layers than the sur¬ 
rounding locations. Synthesized brightness at mm wavelengths 
does not necessarily follow the effective formation height pro¬ 
file, but there is a tendency for local brightenings above the field 
concentrations (see curves of the same type in Fig. 5a-c and 5d). 
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Fig. 7. Spatial location (x,y)=(-8.1,0.9) Mm, marked ”1” in Fig. 6. (a) Electron temperature (solid curve), electron number density (dashed curve), 
and magnetic field strength (dotted) as a function of height, (b) Contribution functions to emerging intensity (normalized to their maximum) as 
labeled in the figure. The horizontal line marks 10% of the highest value for comparison with Fig. 6. (c) Resulting brightness spectrum. Stars mark 
the four wavelengths analyzed in detail. 




Fig. 8. Same as in Fig. 7 for the spatial location (-5.2,0.9) Mm, marked ”2” in Fig. 6. 




Fig. 9. Same as in Fig. 7 for the spatial location (0.7,0.9) Mm, marked ”3” in Fig. 6. 


In Fig. 6 we show the range of geometrical heights that con¬ 
tribute to the emerging intensity as a function of location along 
the vertical slice from Fig. 5. The heights marked in Fig. 6 
with blue contours at a level of 10% of the CF maximum can 
serve as a measure of the height range sampled by each of the 
four analyzed wavelengths. For comparison, we plot the corre¬ 
sponding effective formation height (red solid line) at each wave¬ 
length. Owing to the complex thermal structure, a broad height 
range is involved in emitting radiation at all considered wave¬ 
lengths, resulting in quite complicated contribution functions 
(see Sect. 4.4 for details). The figure demonstrates the limita¬ 
tions of using effective heights to assign the synthesized bright¬ 
ness to the thermal structure of the atmosphere. At a number of 


locations the contribution achieves significant values over two 
distinct height ranges with double-peaked contribution functions 
(see, e.g., A= 1 mm in Fig. 8 and T=3.6 mm in Figs. 9 and 10). 

4.4. Analysis of individual profiles 

In this section, we analyze the formation of the mm brightness 
spectra at five locations along the XZ-slice that represent differ¬ 
ent features in the solar atmosphere and display a broad range 
of behavior. In Fig. 5 these locations are marked by dashed ver¬ 
tical lines, and in the top panel of Fig. 6 they are additionally 
labeled with numbers from 1 to 5. The properties of the spa¬ 
tially resolved atmospheres (electron temperature, electron num- 
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Fig. 10. Same as in Fig. 7 for the spatial location (9.5,0.9) Mm, marked ”4” in Fig. 6. 




Fig. 11. Same as in Fig. 7 for the spatial location (10.2,0.9) Mm, marked ”5” in Fig. 6. 


ber density, magnetic field strength) as a function of geometri¬ 
cal height, contribution functions at four wavelengths and corre¬ 
sponding brightness spectra are shown in Figs. 7, 8, 9, 10, and 
11, for features 1 to 5. Note that for better comparison the scales 
remain unchanged from one figure to the next. 

- The spatial location (x, y) = (-8.1,0.9) Mm (marked "1" in 
Fig. 6) is detailed in Fig. 7: At the solar surface the location 
corresponds to an intergranular feature with a rather weak 
magnetic field (< 100 G). Its temperature distribution is char¬ 
acterized by a low-lying transition region (at 1500 km). The 
electron density manifests a localized increase (by an order 
of magnitude) at around 1400 km. The contribution functions 
at 0.4 mm and 1 mm are very sharply defined, and the con¬ 
tributing formation heights are in the range 600 - 700 km, 
where temperature decreases with height, leading to a tem¬ 
perature minimum at around 800 km. At 4=3.6 mm the con¬ 
tribution function is double peaked, with the main peak at 
lower height covering the extended temperature minimum 
region, and the narrow, weak secondary peak in the chro¬ 
mosphere and lower transition region. As a consequence, the 
synthesized brightnesses are extremely low (around 4000 K) 
and brightness temperatures at 1 mm and 3.6 mm are lower 
than that at 0.4 mm. Radiation at 10 mm is formed over a 
range of heights with a sharp and clear maximum contribu¬ 
tion from around 1450 km. The resulting brightnesses are 
also significantly lower than the average values for this wave¬ 
length. 

- Spatial location (x, y) = (-5.2,0.9) Mm (Fig. 8, marked ”2” 
in Fig. 6): Here B reaches 1900 G in the photospheric in¬ 
tergranular lane. At short wavelengths the contribution func¬ 
tions are double peaked, at 0.4 mm the main peak is at the 
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height of the temperature minimum with a secondary peak 
in the low chromosphere, whereas at 1 mm the main peak 
moves to chromospheric heights, but a small contribution 
from the temperature minimum region remains. At longer 
wavelengths the contribution heights move higher into the 
chromosphere, the contribution functions are sharply defined 
with a single peak close to the boundary between the chro¬ 
mosphere and the TR. The resulting mm brightnesses are sig¬ 
nificantly higher than the average values given in Table 1 and 
increase almost linearly with wavelength. 

- Spatial location (v, y) - (0.7,0.9) Mm (Fig. 9, marked ”3” 
in Fig. 6): This location corresponds to the central part of 
the loop connecting the two main field concentrations and is 
characterized by complicated temperature and density pro¬ 
files showing a number of local enhancements along the line 
of sight. This leads to a compound contribution function at 
3.6 mm. The brightness spectra achieve average values. 

- Spatial location (v, y) = (9.5,0.9) Mm (Fig. 10, marked ”4” 
in Fig. 6): This is a peculiar location at the boundary be¬ 
tween an intergranular lane and a granule. It is characterized 
by a bubble of very hot (basically coronal temperature), low- 
density gas at a height of about 2000 km, somewhat like a 
piece of corona, separated from the main corona. This fea¬ 
ture, due to low electron density, produces low opacity and 
does not contribute to the emerging intensity. As a result, the 
contribution function at 3.6 mm has peaks on either side of 
this abrupt temperature enhancement. Radiation at 0.4 mm 
and 1 mm forms below these layers and is not affected by 
this inhomogeneity, while radiation at 10 mm is formed only 
above it. The brightness at 0.4 mm and 1 mm is typical for 
these wavelengths. The 7^ at 3.6 mm and 10 mm exceeds the 
average values. 
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- Spatial location (x,y) = (10.2,0.9) Mm (Fig. 11, marked 
”5” in Fig. 6): This is an example of the atmosphere above a 
granule. A steep transition region at 2000 km height results 
in a narrow contribution function for 10 mm and in an ex¬ 
tremely low synthesized brightness. The radiation at shorter 
wavelengths receives its main contribution from the heights 
where the temperature is below 5000 K, and as expected, it 
also shows very low brightness. Both this location and the 
previous one display signs of a magnetic canopy with a base 
at around 500-700 km. 


0.4 mm 1 mm 



3 mm 4.5 mm 




Tgas (kK) Tgas (kK) 


Fig. 12. Density scatter plots (2D histograms) of the brightness tem¬ 
perature T b versus gas temperature T gas taken at the effective formation 
heights for 0.4 mm, 1 mm, 3 mm, and 4.5 mm. Darker color indicates 
more pixels in the bin. Solid lines denote T b = T gas . Blue contour (dark 
gray in the printed version) encloses 95% of the spatial locations with 
weak magnetic field (|5| < 30 G at the effective formation height), red 
contour (light gray in the printed version) includes 95% of the locations 
of strong magnetic field (|5| > 30 G at the effective formation height). 



Fig. 13. Brightness temperature at 3.6 mm (solid), electron tempera¬ 
ture at the effective formation height (dashed) along the X-profile at 
Y=0.9 Mm. The vertical dashed lines are the same as in Figs. 5 and 6. 


4.5. Thermal diagnostics 

In this section we investigate the potential of submm/mm bright¬ 
ness spectra to diagnose chromospheric plasma. We present cor¬ 
relations between brightness temperature and the properties of 
the atmosphere, in particular its thermal structure. 


In Fig. 12 we correlate the synthetic mm brightness with the 
model temperature taken at the heights corresponding to the ef¬ 
fective formation heights at wavelengths 0.4 mm, 1 mm, 3 mm, 
and 4.5 mm. Because of the many pixels in the snapshot, we 
show the correlations by binning the data into 2D histograms. In 
the density scatter plots we distinguish between a weak magnetic 
field with \B\ < 30 G at the effective formation height (darker 
gray contours in the printed version and blue contours in the on¬ 
line version), and an enhanced magnetic field with \B\ > 30 G 
(light gray contours in the printed version, red contours in the 
online version). The scatter seen in Fig. 12 illustrates the influ¬ 
ence of the extended formation height ranges, which instead of 
sampling the temperature at a fixed geometrical height, effec¬ 
tively mixes contributions from different layers. 

For regions with large B , the largest differences between T gas 
and T b occur at the high and low end of the temperature range 
sampled at each A. In the cool part of the range mm brightness 
tends to overestimate T gas , while when T gas is high the synthetic 
brightness underestimates it. For weak magnetic features a wider 
’’scatter rectangle” around the T gas = T b curve in the scatter plots 
is seen. When confining the analysis to fields with \B\ > 30 G, 
the scatter becomes much smaller and mm brightness reproduces 
temperature at different chromospheric heights relatively well. 
The agreement for \B\ > 30 G is better for longer wavelengths, 
which has two reasons: applying a 30 G magnetic field strength 
threshold at the formation height means that we consider fewer 
regions at longer A , since on average B decreases with height. In 
the more strongly magnetized parts of the atmosphere contribu¬ 
tion functions at longer wavelength show a simple form with one 
peak, see Fig. 8, so that the effective formation height provides a 
good measure of the real formation height. 

The relation between brightness temperature and atmo¬ 
spheric temperature can also be seen in Fig. 13, where we 
show brightness at 3.6 mm (solid), the longest of the cur¬ 
rently available ALMA wavelengths, as a function of position 
along the analyzed XZ-cut, together with the model tempera¬ 
ture (dashed) taken at the heights corresponding to the effective 
formation heights. In the T gas profile there are excursions to val¬ 
ues lower and higher than the synthesized T b . Furthermore, at 
X = +9.5 Mm, there is a huge discrepancy between the two 
quantities. It corresponds to an intergranule-granule boundary. 
This location was analyzed in Fig. 10. At this location the effec¬ 
tive height evaluated from the centroid of the contribution func¬ 
tion does not reflect the real formation heights because of the 
complex form of the CF that has two peaks of similar strength 
with negligible contribution from in between. For strong mag¬ 
netic elements (e.g., spatial location X = -5.2 Mm) the agree¬ 
ment between synthetic brightness and atmospheric temperature 
is nearly perfect. 

4.6. Effect of spatial smearing 

Depending on array configuration (maximum baseline) and 
wavelength, ALMA can achieve a spatial resolution in the range 
0.005" - 5". In Table 2 we list the estimates of the FWHM of the 
synthesized beam (point spread function), which is the inverse 
Fourier transform of a (weighted) u-v sampling distribution. To 
account for this finite resolution, we spatially smeared the syn¬ 
thetic brightness maps by convolving with a Gaussian kernel of 
corresponding FWHM to mimic the instrumental profile. Note 
that we did not change the pixel size, so that in the following 
plots we oversample the spatially smeared data sets. 

The effect of spatial smearing on the mm brightness is il¬ 
lustrated in Fig. 14 for three wavelengths: 0.1 mm, 1 mm, and 
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Fig. 14. Spatially smeared synthetic brightness at 0.4 mm (top row), 1 mm (middle row), and 3 mm (bottom row) at a spatial resolution of 0.2", 
0.4", 1" and 4". The same color scale is used for all panels in a row. The field size is 24 Mm x 24 Mm. 



brightness (K) 


1 mm 



brightness (K) 


3 mm 



brightness (K) 


3.6mm 0.2" 


0.4" 




Tgas (kK) 

1 " 


Tgas (kK) 
4" 




Tgas (kK) 


Tgas (kK) 


Fig. 16. Density scatter plots of brightness temperature at 3.6 mm at 
resolutions of 0.2", 0.4", 1" , and 4" versus gas temperature T gas taken 
at the effective formation heights. Darker color indicates more pixels in 
the bin. The data are restricted to the locations of strong magnetic field 
(|5| > 30 G) at the effective formation height. Red contours enclose 
95% of the selected pixels. 


Fig. 15. Spatial resolution effect on the brightness histogram. Top: 
brightness histograms at 0.4 mm for a spatial resolution of 0.2" (solid), 
0.4" (dotted), 1" (dashed), and 4" (dot-dashed). Middle: the same for 
1 mm. Bottom: the same for 3 mm. 


3 mm, and for four values of FWHM: 0.2", 0.4", 1" , and 
4". While at 0.2" much of the fine structure of the original- 
resolution images is preserved, it is increasingly lost as the res¬ 
olution is reduced further. At the same time, the intensity con¬ 
trast is rapidly reduced, as can be deduced from the intensity 
histograms in Fig. 15 plotted for the same three wavelengths 
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Table 2. Spatial resolution of the ALMA interferometer as a function 
of the longest array baseline for a number of wavelengths. 


A, mm 

FWHM(") 
maxbase= 150 m 

FWHM(") 
maxbase= 16 km 

0.4 

0.5 

0.005 

1.0 

1.3 

0.01 

3.0 

4 

0.04 

4.5 

6 

0.06 


and from Table 3. In Table 3 we list brightness variations T^ ms 

jrms 

and relative brightness temperature contrast as a function 
of ALMA resolution for comparison with the original values in 
Table 1. 

The difference between original brightness and brightness at 
a resolution of 0.2" is negligible. This is expected because of 
the rather large grid spacing of 48 km of the original simula¬ 
tion, which makes the spatial resolution of the original simu¬ 
lation similar to the degraded one. At a resolution of 0.4" the 
finest structure starts to be washed out (but insignificantly), and 
the most extreme brightness values are no longer present in the 
histogram (see dotted curves in Fig. 15). At all wavelengths the 
relative brightness contrast changes only very little from its orig¬ 
inal values. However, there is a significant change at 1" res¬ 
olution: the contrast begins to be reduced significantly (it de¬ 
creases to 70% of original contrast at the shortest wavelengths), 
which is visible as a narrowing of the intensity histogram (see 
dashed curves in Fig. 15), the fine structure is already smeared 
out, while the larger scale pattern is well distinguishable and 
closely reflects that in the original image. Finally, at 4" reso¬ 
lution, which corresponds to the best resolution of the solar in¬ 
terferometric observations with CARMA (Combined Array for 
Research in Millimeter-wave Astronomy) currently available, all 
the fine structure is lost and the contrast has dropped dramati¬ 
cally. 

The good agreement between synthetic mm brightness and 
gas temperature evaluated at the effective formation height 
steadily degrades as the spatial resolution is decreased. In Fig 16 
we show in the form of a 2D histogram how it develops for 
4=3.6 mm from the original dependence to a spatial resolution 
of 4" for the computed T\ maps. Note that the electron temper¬ 
ature at the effective formation height is not spatially degraded, 
since we wish to determine how much of the true local temper¬ 
ature can be deduced from degraded data. The range covered by 
Tb decreases with decreasing resolution, the electron tempera¬ 
ture does not, so that the relationship starts to deviate from the 
ideal line, mainly at low temperatures. Nonetheless, the diagnos¬ 
tics remains fairly robust up to 1" resolution, particularly for the 
somewhat higher temperatures. At 4" resolution, however, the 
observed Tb carries almost no information on the original elec¬ 
tron temperature in the fine-scale atmosphere. 


5. Discussion and conclusions 

Using realistic simulations of the quiet Sun, we studied sim¬ 
ulated brightness at millimeter and submillimeter wavelengths 
and how it is affected by the instrumental resolution of ALMA. 

We conclude that millimeter brightness provides a reason¬ 
able measure of the chromospheric thermal structure at the 
height at which the radiation is formed. We found that the bright¬ 
ness temperature is a reasonable measure of the gas temperature 


at the effective formation height, defined as the height corre¬ 
sponding to the centroid of the contribution function, in spite 
of considerable scatter. The scatter is significantly reduced when 
very weak magnetic fields are avoided. Our results indicate that 
although instrumental smearing reduces the correlation between 
brightness and temperature, Tb can still be used to diagnose the 
electron temperature up to a resolution of 1". Another effect of 
finite spatial resolution is that it smooths out the more extreme 
values. 

When comparing images, spatial resolution emerges as a 
critical problem for chromospheric observations with ALMA. 
It will be important to achieve the highest possible resolu¬ 
tion because the fine details are increasingly lost as the reso¬ 
lution decreases. We cannot conclude from the current simula¬ 
tion whether 0.1" is a sufficiently high resolution because of the 
rather large grid size of 48 km=0.064" of these simulations. Sim¬ 
ulations with a finer grid may reveal more fine structure. 

The often complex contribution functions of the mm radi¬ 
ation imply that the brightness temperatures represent the in¬ 
tegrated physical state of an extended height range, sometimes 
also over two distinct, well-separated height ranges. However, 
multiwavelength observations with ALMA with a narrow step 
in wavelength can provide chromospheric images formed over 
subtly different height ranges. This will make tomography of the 
chromospheric thermal structure feasible (called ’’volume imag¬ 
ing” by Wedemeyer-Bohm et al. (2007)). This conclusion is 
confirmed by Fig. 17, which displays density scatter plots of 
brightness temperatures at three wavelengths corresponding to 
different ALMA frequency bands. Neighboring bands (Figs. 17a 
and 17c) show less deviation from diagonal lines, depicting ex¬ 
pected identical brightness for the two wavelengths and thus 
carry some overlapping information, in contrast to the bands far¬ 
ther away (Fig. 17b), which show larger scatter. 

ALMA will provide a unique opportunity to probe the solar 
atmosphere from the height of the classical temperature mini¬ 
mum through the chromosphere to just below the transition re¬ 
gion (see Fig. 5). The shortest currently available wavelengths 
at ALMA will examine the temperature minimum region in the 
solar atmosphere, while in the range from 1 mm to 3 mm we can 
gain access to the middle chromospheric heights (at 1000 km- 
2000 km). At the longest wavelengths, planned for the future 
development of ALMA, the transition region becomes acces¬ 
sible. Moreover, solar chromospheric science can profit from 
the ALMA observations at submillimeter and millimeter wave¬ 
lengths used in combination with other chromospheric diagnos¬ 
tic such as Ca n K, Mg n k (Leenaarts et al. 2013a,b; Danilovic 
et al. 2014; Riethmuller et al. 2013), Ca n 854.2 nm (de la Cruz 
Rodriguez et al. 2013), or Ha (Leenaarts et al. 2012). 

As the next steps of this study, we are planning to study diag¬ 
nostics of magnetic fields at the heights where submm-mm radi¬ 
ation is formed. We will discuss simulated circular polarization 
in the context of successful solar ALMA polarization measure¬ 
ments. We will reconsider the assumptions made in this work 
for computing the submm-mm radiation, in particular the use of 
quasilongitudinal approximation. Next we plan to go beyond a 
single snapshot to investigate the time dependence and the nec¬ 
essary time resolution of ALMA observations. 

In this paper we used one of the most realistic simulations of 
the solar chromosphere currently available. However, as reported 
by Leenaarts et al. (2013b), there are two processes of particular 
importance for the chromospheric thermal structure that are not 
included in the numerical simulations analyzed here. These refer 
to the non-equilibrium ionization of helium and the effect of par¬ 
tial ionization on chromospheric heating by magnetic fields (see 
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j^rms 

Table 3. RMS variation T™ ls and relative brightness temperature contrast for a number of analyzed wavelengths as a function of the ALMA 
spatial resolution 


Resolution 
(arc sec) 

0.4 mm 

r rrms 

1 b 

1.0 mm 

(K) 

3 mm 

4.5 mm 

0.4 mm 

'T'rms 

1 b 

<T b > 

1.0 mm 3 mm 

4.5 mm 

0.2 

496 

811 

1311 

1465 

0.11 

0.17 

0.22 

0.22 

0.4 

442 

738 

1220 

1369 

0.10 

0.15 

0.20 

0.20 

1.0 

339 

591 

1041 

1192 

0.08 

0.12 

0.17 

0.18 

4.0 

153 

310 

667 

819 

0.03 

0.06 

0.11 

0.12 


1 mm vs. 0.4mm 



Tb(10 3 K) 0.4mm 


3.6mm vs. 0.4mm 



Tb(10 3 K) 0.4mm 


3.6mm vs. 1mm 



Tb(10 3 K) 1mm 


Fig. 17. Density scatter plots of brightness temperatures T b for T=0.4 mm, 1 mm, and 3.6 mm. Darker color indicates more pixels in the bin. Blue 
contours enclose 95% of the pixels. Solid lines denote T= T* 2 . 


Leenaarts et al. 2013b, and references therein). Studies based on 
simulations of regions with different amounts of magnetic flux 
would also be useful, as well as simulations with a smaller grid 
size. 
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